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Abstract. The potential of the Surface Brightness Fluctuations (SBF) method to determine the membership 
of dwarf elliptical galaxies (dEs) in nearby galaxy clusters is investigated. Extensive simulations for SBF mea- 
surements on dEs in the I band for various combinations of distance modulus, seeing and integration time are 
presented, based on average VLT FORS1 and FORS2 zero points. These show that for distances up to 20 Mpc 
(Fornax or Virgo cluster distance), reliable membership determination of dEs can be obtained down to very faint 
magnitudes — 10 < M\ < — 12 mag (/^o(V) ~ 25 mag arcsec -2 ) within integration times of the order of 1 hour 
and with good seeing. Comparing the limiting magnitudes of the method for the different simulated observing 
' conditions we derive some simple rules to calculate integration time and seeing needed to reach a determined 

, limiting magnitude at a given distance modulus for observing conditions different to the ones adopted in the 

■ simulations. Our simulations show a small offset of the order of 0.15 mag towards measuring too faint SBF. It 
is shown that this is due to loss of fluctuation signal when recovering pixel-to-pixel fluctuations from a seeing 
convolved image. To check whether our simulations represent well the behaviour of real data, SBF measurements 
for a real and simulated sample of bright Centaurus Cluster dEs are presented. They show that our simulations 
are in good agreement with the achievable S/N of SBF measurements on real galaxies. 

|' Key words, galaxies: clusters - galaxies: dwarf - galaxies: fundamental parameters - galaxies: luminosity function 

^ \ - galaxies: distances and redshift - techniques: photometric 

■ 1. Introduction a is very useful for testing models of galaxy formation. 



There are two fundamental steps involved in determining 
1.1. The faint end of the galaxy luminosity function the faint end of $(M) in galaxy clusters: first; finding the 
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Dwarf elliptical galaxies (dEs) are the most numerous dwarf S alax y candidates; second, verifying that they are 

03 type of galaxies in the nearby universe, especially in clus- cluster members. 

ters. This statement has been well established since the For the Local Grou P> $ ( M ) has been determined down 

advent of CCD detectors and the building of telescopes to M v - " 9 ma § ( e -§- Mateo QUI Pritchet MM Van 

with 4-8m diameter that enabled observers to detect low den Ber S h 2000 >' suggesting a ~ -1.1. Local Group 

surface brightness (LSB) objects substantially fainter dwarf galaxies can nowadays readily be resolved into 

than the night sky. With the improvement of observing sin § le stars with HST and / or active °P tics techniques, 

facilities, the emphasis has over the last decade switched Thus ' their Stance can be determined; the second step 

from detecting faint dwarf galaxies to quantifying well in establishing $(M) is quite easy to perform. The first 

their properties and frequencies. Most of the times dEs ste P> nndin g them ' is more difficult due to their large 

are investigated in galaxy clusters, because for a cluster angular extent and small contrast against the stars of 

the distance and therefore the approximate angular size the Milk y Wa y- The latest discoveries of more and more 

of candidate dwarf galaxies is known. ^ dSphs (e.g. Armandroff et al. QUI Whiting et al. 

One of the most important statistical tools in investigat- HS23 raise the question of how complete the Local Group 

ing galaxy populations is the galaxy luminosity function sample is. 

$(M), describing the frequency of galaxies per magnitude The opposite is true for nearby galaxy clusters. Here, 
interval. The knowledge of its logarithmic faint end slope nndin g candidate dwarf spheroidals is quite straightfor- 
ward when performing deep enough photometry, but 
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example, the brightest red giants of an early type galaxy 
at the Fornax cluster distance (19 Mpc, Ferrarese et al. 
HUUty have V ~ 29.4 mag (Bellazzini et al. I3UUTJ) . The 
confirmation of candidate dwarf spheroidals as cluster 
members must consequently be based on morphology and 
is therefore subject to possible confusion with background 
LSB galaxies. One depends on statistical subtraction 
of background number counts to estimate the faint end 
slope. Due to the generally low number counts, the 
Poisson error involved in this statistical subtraction 
constitutes a major source of uncertainty in determining 
a, especially in magnitude-surface brightness bins where 
the contribution of background galaxies is of the order of, 
or higher than, that from cluster galaxies. 
The majority of studies that have determined $(AI) in 
nearby galaxy clusters (e.g. Sandage et al. 1198 5. Trentham 
et al. I2UUT1 I^OO^al 12002 bl and Hilker et al. |2HU3|) suggest 
a logarithmic faint end slope of —1.0 < a < —1.5, being 
in substantial disagreement with CDM theory (Press 
& Schechter I1974jl . which predicts a ~ -2.0 for the 
initial galaxy luminosity function (Kauffman et al. 2000 ) . 
Other authors, such as Phillips et al. 1|1998H for the Virgo 
cluster and Kambas et al. I|2000|l for the Fornax cluster, 
suggest a significantly steeper faint end slope of a ~ — 2. 
This large discrepancy shows that much care must be 
taken when assigning cluster membership to galaxies for 
which no direct distance measurement is available. The 
poisson statistics involved, especially when subtracting 
background number counts, can lead to different authors 
obtaining very different results for the same cluster. 
Various methods for distance determination that can 
be applied to brighter galaxies outside the Local Group 
are not suited for the faintest dEs: standard candles 
such as SN la or Cepheids are very rare; radial velocity 
measurements, if possible, are very time consuming due 
to the large fields that have to be covered. 



1.2. Distances to galaxies with the SBF Method 

An intriguing possibility to unambigously determine 
cluster membership of large samples of dEs in nearby 
clusters is deep wide field imaging and application of 
the surface brightness fluctuations (SBF) method. The 
SBF-method was first described by Tonry & Schneider 
(1988). SBF are caused by the fact that on a galaxy 
image the number of stars in each seeing disc is finite 
and therefore subject to statistical fluctuations. The 
amplitude of these fluctuations relative to the underlying 
mean surface brightness is inversely proportional to the 
number of stars. As the number of stars per unit angle 
increases quadratically with distance, the amplitude of 
the SBF is inversely proportional to distance and can 
therefore serve as a distance indicator. 
As it is well known, the integrated light of an entire 
stellar population is dominated by the light emitted 
from the giants. In the case of an old population, only 



the red giants contribute. In the case of a younger 
population, there is also a contribution from the blue 
super giants. In the context of the SBF-method one can 
treat the whole stellar population as consisting only of 
stars having the mean luminosity weighted luminosity 

— f n(L)LdL 

L = <*-* of this population. The absolute surface 

J n(L)dL t * 

brightness fluctuation magnitude M is then defined as 
M = —2.5 x log(L) + c with c being the zero point of 
the photometric system. To obtain the distance modulus 
of a galaxy with the SBF-method, one directly measures 
the apparent SBF magnitude m and derives M from 
a distance- independent observable, usually a colour. 
Measuring the apparent magnitude and deriving its 
absolute value from a distance-independent observable 
is a procedure common to many distance determination 
methods, e.g. the P-L relation for Cepheids. 
Tonry et al. 1TTW71 I2UUU1 EMTjl have carried out an 
extensive survey to measure SBF for bright elliptical and 
spiral galaxies in 22 nearby galaxy groups and clusters 
within 40 Mpc. Using distances derived from cepheids or 
SN la in the respective galaxies, they obtain a relation 
between the absolute fluctuation magnitude Mi and the 
dereddened colour (V — I), determined in the colour 
range 1.0 < (V - I) < 1.3: 

Mi = -1.74 + 4.5 x {{V -I)- 1.15) mag (1) 

Up to now the SBF-method has only been applied to 
small samples of nearby dEs (e.g. Bothun 119911 Jerjen et 

ai. nnsHi Enni and onniii • 



1.3. Aim of this paper 

In this paper, we focus on the potential of the SBF 
method to determine cluster membership of dEs in 
nearby clusters. We measure SBF of simulated sets of dEs 
at 3 different distance moduli between 7.5 and 43 Mpc 
for various combinations of observing times and seeing. 
The paper is structured as follows: In Sect.El the impact 
of the age-metallicity degeneracy on SBF magnitudes is 
discussed. In Sect. 01 the technical details of simulating 
and measuring SBF are described and the properties of 
the different sets of simulated dEs are shown. In Sect. 
the results of the SBF measurements for all simulated sets 
are presented and the limiting absolute magnitudes for 
cluster membership determination is discussed. A validity 
check of the simulations is presented by comparing real 
and simulated SBF data. We finish this paper with the 
conclusions in Sect. |5j 



2. Deriving Mi from (V - I) 

In the context of the SBF-method, the distance modulus 
of a galaxy is given by the difference between apparent 
and absolute fluctuation magnitude (toi — Mi). To 
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estimate the reliability of the method, one must know 
both the accuracy in measuring m\ at the cluster distance 
and the uncertainty in deriving M j for a dE with a given 
(V-I). _ 

In Fig.^ theoretical values for Mi are plotted vs. (V — /) 




0.6 0.8 1 



(V-I) [mag] 

Fig. 1. Theoretical absolute fluctuation magnitude M\ 
plotted vs. the colour (V — I) for a set of different stellar 
populations, as taken from Worthey (1994). The larger the 
symbols, the older the stellar population. Largest symbols: 
17 Gyr; intermediate size symbols: 12 Gyr; small symbols: 
8 Gyr (except for the open triangles, see below). Solid 
circles: Single burst populations with Fe/H between -1.7 
(-1.3 for 17 Gyr) and dex (left to right). Open circles and 
asterisks: 10% and 30% mixing of the single burst popu- 
lations from above with a 5 Gyr population of Fe/H=0. 
Open triangles: Single burst population of 1.5 Gyr and 
Fe/H > -1.3. The solid line is defined by Relation 

for a set of old to intermediate age stellar populations 
with a wide range of metallicities, taken from Worthey 1 
(1994). Equation is indicated as well. Note that Tonry 
et al. determined equation only for 1.0 < (V — /) < 1.3 
as they used brighter and therefore redder galaxies. The 
theoretical and observed values for M\ agree quite well, 
although for (V — I) > 1.15 the slope of equation is 
slightly shallower than for the theoretical values. 
Fig. ^ shows that for (V — I) > 1.0, Mj depends mainly 
on (V — I) and very little on the age-metallicity combi- 
nation of the underlying stellar population. However, for 
(V — I) < 1.0 the spread in M\ at one given (V — I), i.e. 
the effect of the age-metallicity degeneracy, rises signif- 
icantly. At (V-I) ~ 0.85 this leads to an uncertainty 
of up to 0.5 mag in relating M\ to (V — I). Jerjen et al. 

1 http://astro.wsu.edu/worthey/dial/ dial_a_pad.html 



(1998, 2000) have found an analogous spread of the order 
of 0.5 mag from i?-band SBF measurements for blue 
dEs in the Sculptor and Centaurus A group. Additional 
colour information can help to reduce this uncertainty; 
ej^, at constant [V - I) = 0.9, d~Mi/d(V - K ) ~ 10 and 
d~Mi/d(U - V) ~ 5 (Worthey 1994). Precise photometry 
in U, V and I would therefore allow a precision of 0.1 to 
0.2 mag in deriving M\ from theoretical models. 
While the age-metallicity degeneracy is an important 
source of uncertainty for distance measurements to 
blue field galaxies, it becomes useful as a relative 
age-metallicity indicator for blue cluster galaxies, as 
most of the time the cluster is separated from back- 
ground/foreground galaxies by significantly more than 
0.5 mag in distance modulus. Spectroscopic surveys of 
the Fornax cluster (Drinkwater et al. 2000, Hilker et al. 
1999) revealed a significant gap in radial velocity between 
Fornax members and background galaxies corresponding 
to ~ 3 mag in (m — M) . Confusion of a blue and young 
background galaxy with a blue and old cluster galaxy is 
therefore very unlikely to happen. 

Note that dlrr candidate members, to which the SBF- 
method is difficult to apply due to their irregular shape, 
can be distinguished morphologically in a straightforward 
manner from blue young background galaxies like anemic 
spirals. For clusters with a significant fraction of dlrrs, 
this morphological cluster membership assignment can 
complement the assignment based on SBF-distances for 
the smoothly shaped dE candidates and allow derivation 
of $(M) for the entire dwarf galaxy population. 



3. Simulating and measuring Surface Brightness 
Fluctuations in dEs 

We have simulated sets of dEs in the /-band with 3 
distance moduli 29.4, 31.4 and 33.4 mag, corresponding 
to 7.6, 19 and 48 Mpc. This range was chosen to include 
distances to the more nearby groups like Leo I (10 Mpc) 
as well as to the more distant clusters like Centaurus 
and Hydra (33 to 33.5 mag). Note that 31.4 mag is the 
approximate distance modulus to Fornax and Virgo. The 
integration time was 3600 seconds, the gain was 1 and 
the zero point 27.0 mag. The latter value is a mean of 
the VLT FORS1 and FORS2 zeropoint for imaging in 
the /-band when including an averaged colour term and 
extinction coefficient. The pixel scale was 0.2"/pixel, the 
image size 2048x2048 pixel. For each distance modulus, a 
set with 0.5" and 1.0" seeing was simulated. Additionally, 
for 31.4 mag distance modulus and 0.5" seeing, four 
different integration times were adopted, namely 900, 
1800, 3600 and 7200 seconds. Note that varying the 
integration time from t\ to ti can also be considered as 
keeping the integration time fixed and adding 2.5 log -jj^ 
to the zero point. 

The photometric properties of our simulated dEs are 
derived from Hilker et al.'s (2003) values found for Fornax 
cluster dEs. For a given absolute magnitude My, the 
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colour-magnitude- and surface brightness-magnitude 
relation from Hilker et al. (2003) is used to obtain (V — I) 
and /x (V). For example, a dE with My = — 10 mag will 
have (V — I) = 0.86 mag and Ho(V) = 25.9 mag arcsec -2 
while a dE with M v = -15 mag will have (V - I) = 1.04 
mag and /Lto(V) = 22.7 mag arcsec -2 . An exponential 
intensity profile of the form I(r) — Iq x exp(~r/ro) was 
adopted, with ro calculated from My, ^to(V) and the 
adopted distance modulus. The ellipticity was chosen as 
zero. Globular Cluster (GC) systems are included, with a 
specific frequency of Sn = 5 (Miller et al. I1998[) for all 
galaxies, an absolute turnover magnitude Mj — — 8.5 mag 
(Kundu et al. 12001(1 and the projected spatial GC density 
following the galaxy light distribution. 
The final output of a SBF measurement procedure is 
the apparent surface brightness fluctuation magnitude 
m (in our case mi). This is equivalent to the luminosity 
weighted average apparent luminosity of the observed 
stellar population. As mi = Mi + (to — M) , with (to — M) 
being the adopted distance modulus of the simulated dE, 
Mi determines the SBF amplitude at a given distance 
modulus. Mi was adopted as a function of (V — I) 
according to Tonry's equation (JIJ for (V — I) > 1.0. For 
(V — I) < 1.0, it was decided to split the sample into two 
halves, since it is not known whether Tonry's equation 
also holds for galaxies bluewards of this limit. For 50% 
of our galaxies, Mi was calculated according to equa- 
tion Q). For the other 50% of our galaxies, Mi was kept 
constant at Mi((V — I ) = 1.0)) = -2.4 mag. Doing so 
we acknowledge that there is a significant age-metallicity 
degeneracy for (V — I) < 1.0 (see Worthey's models). 
For the bluest galaxies simulated, with (V — I) ~ 0.8 
(My — —8.5 mag), this implies a range in Mi of about 
0.9 mag between the two simulated samples. The effect 
for the simulations is that, on average, the SBF signal is 
weaker than if all dEs with (V — I) < 1.0 were simulated 
according to equation 0J. 

To allow for varying seeing and integration times, the 
background field had to be created artificially. In a 
real background field obtained with VLT FORS1 in the 
/-band at 3000 sec integration, we fitted a power law 
distribution of the form n(m) = A x lQ"ii m ^ m o) to 
the magnitude distribution of the objects detected by 
SExtractor down to the completeness limit of I ~ 25 mag. 
The fitted values were A = 14900/degree 2 , 7 = 0.305, 
mo = 22. According to this distribution, an artificial 
object field was created with the IRAF task mkobjects 
in the ARTDATA package, with 7=27 mag as the faint 
limiting magnitude. Seeing and integration time were 
chosen as needed for the simulations. The sky brightness 
was adopted as 19.9 mag arcsec -2 in I, which holds 
within 3 days before and after new moon. To simulate 
large-scale flat-field effects, the background fields were 
multiplied by the normalized sky-map obtained from 
applying SExtractor to a VLT FORS1 flat field image. 
Into each 2048x2048 pixel field, 16 dEs were imple- 
mented. For each set with constant seeing, distance 
modulus and integration time, six fields were created, 



each with a different (random) spatial distribution of 
the background objects. This means that 96 dEs were 
simulated for each set. 



3.1. Simulation of Surface Brightness Fluctuations 

For each modelled pixel with a given distance r to the 
galaxy center, first the number of stars corresponding to 
the surface brightness fj,(r) of the exponential profile was 
calculated: 

AWs/pixelW = IQ-^xiMrHm-M)-^)) x p 2 (2) 

with p being the pixel scale. Then, a random number 
N* was chosen within a Poisson distribution centered on 
-^stars/pixel ( r )- The intensity adopted at that pixel was 
then defined as 

AT* 

I(pixel,r) = W -0AX(^r)-ZP) x p 2 x f* (3) 

-^stars/pixel v ) 

with ZP being the zero point, 27 mag in our case. 
The implementation of the SBF is apparently achieved 
by multiplying by — — — -j-pj. This means that along 

an isophote with radius r, the intensity I(r) has a 
pixel-to-pixel rms of . I ^ The image with the 

implemented SBF was then convolved with a Moffat 
seeing profile, which was modelled out to 7 times the 
FWHM. Finally, Poisson noise with rms=v / 7 was imple- 
mented. Once modelled, the galaxies were added onto the 
artificial background fields. 



3.2. Measurement of SBF 

To measure the SBF of a simulated dE, the following steps 
were undertaken: 

1. Create object map with SExtractor of the whole image 
(containing 16 dEs) 

2. Mask the dEs on the object map, subtract this image 
from original image 

3. Create and subtract SExtractor sky map 

4. Determine and subtract local sky level by a curve 
of growth analysis with the IRAF-task ELLIPSE in the 
ISOPHOTE package 

5. Model mean galaxy light with ELLIPSE using a sigma 
clipping algorithm to disregard contaminating sources hid- 
den below the galaxy light, subtract the model 

6. Divide resulting image by square root of the model, cut 
out portion where SBF are measured 

7. Mask out contaminating sources like foreground stars, 
background galaxies and globular clusters. 

8. Calculate the power spectrum (PS) of the cleaned im- 
age 

9. Obtain the azimuthal average of the PS 

10. Fit function of the form 

P{k) = PSF(k) x P + Pi (4) 
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to the azimuthally averaged PS. 

PSF(k) is the PS of the seeing profile, normalized to unity 
at k=0. PSF(k) is determined from a simulated star with 
no close neighbours by fitting a Moffat profile to its PS. 
Pi is the white noise component, proportional to the ra- 
tio between sky and galaxy brightness in the range where 
SBF were measured. It is independent of seeing. Pq is the 
amplitude of the pixel-to-pixel surface brightness fluctua- 
tions, being the zero wavenumber limit of the seeing con- 
volved pixel-to-pixel star count fluctuations, and therefore 
seeing-independent, too. 
It holds that 

mi = -2.5 * /og(P /^exposure) + ZP (5) 

Values at small k (long wavelength) are rejected for the 
fit, as they are often considerably influenced by large-scale 
residuals from imperfect galaxy subtraction and the finite 
width of the image portion used to measure SBF. 
The (seeing independent) S/N of the measurement was de- 
fined as S/N=P /Pi, following Tonry&Schneider ljT3gg|l . 
In the following it will be referred to as canonical S/N. 
Note however that the detectability of SBF decreases 
with increasing seeing: by convolving with the seeing the 
star count pixel-to-pixel fluctuations are smoothed out. 
Additionally, the larger the seeing is, the fewer the in- 
dependent data points per unit angle that are available. 
To take this dependence of the SBF detectability on the 
seeing into account, we defined a modified signal to noise 
S/N* in the following way: 

S/N* = Po/P x y/N^/af (6) 

with iV s d being the number of seeing discs - i.e. indepen- 
dent data points - contained in the image portion where 
SBF are measured, and sf being the smoothing factor 
by which the seeing convolution reduces the pixel-to-pixel 
fluctuations, sf was 5 for 0.5" seeing and 10 for 1.0" see- 
ing, as the effective area of one seeing disc was 25 and 100 
pixels, respectively. 

3.2.1. Examples of simulated dEs 

In Fig. |21 example images and power spectra of simulated 
dEs are shown. The innermost pixels which were neglected 
in the power spectrum fit are especially marked. Note the 
effect of the twice as large seeing of 1.0" in the third ex- 
ample: the width of the seeing power spectrum is about 
half that of the other three examples, where the seeing 
was 0.5". 

Rejecting the innermost pixels is crucial to determine the 
correct SBF amplitude, as low wavenumbers are affected 
by imperfect galaxy subtraction and large scale sky gradi- 
ents in the investigated image. The limiting wavenumber 
beyond which one has to reject pixels has to be determined 
individually for each galaxy, as image dimensions and loci 
and number of contaminating sources change. We have 
adopted the following criterion for deciding which pixels 
to reject or not: if the \ 2 of the fit improves by more than 



a factor of 2 when rejecting the innermost pixel, it is re- 
jected. Then, the same is tested for the second pixel, and 
so on until \ 2 improves by less than a factor of 2. For the 
examples given in Fig. this criterion works fine for the 
lower three power spectra. 

Unfortunately, as illustrated in the upper example, things 
can be more complicated. If only rejecting the inner- 
and outermost data point, the obtained fit fits well to 
wavenumbers smaller than 6, but underestimates the sig- 
nal for wavenumbers between 6 and 10 and overestimates 
the white noise component Pq. When rejecting wavenum- 
bers smaller than 6, the outer part is fit much better. The 
difference in Pi between the two possibilities is consid- 
erable, about 40%. In cases like that, the uncertainty in 
which pixels to reject or not is the major source of error. 
As is noted in Fig. 0s caption, the difference Smi between 
simulated and measured mi obtained when rejecting the 
inner 5 pixels is much smaller than for only rejecting the 
innermost pixel. Therefore, whenever fits to the outer and 
inner part of the power spectrum differed considerably, 
more emphasis was put on fitting well the outer part. 
However, further out in the FT profile the white noise 
component starts to dominate over the PSF FT. As our 
simulations were performed on artificial dwarf galaxies, 
the image portions chosen for the SBF measurements had 
relatively small dimensions of typically 30 to 60 pixel (6 
to 12 "). Therefore the wavenumber range over which the 
amplitude of the PSF is determined is only of the order 
of 10 or fewer independent data points. This small num- 
ber, together with the uncertainty in which wavenumbers 
to disregard or not, is the major source of uncertainty for 
the SBF measurements we performed. For the simulated 
galaxies at 1.0" seeing, the wavenumber range over which 
to perform the fit to the PS is only half of that for 0.5" 
seeing, which means the uncertainty at 1.0" seeing is sig- 
nificantly higher than for 0.5" seeing. 



4. Results of the simulations 

In Figs. 01 to El the results of the simulations are shown. 
Figs. to El show the results for the three different 
distance moduli with 3600 seconds integration time: 
S/N* and Smi are plotted vs. My for the two different 
seeing values. In Fig. the same observables are plotted 
for four different integration times at a fixed 31.4 mag 
distance modulus and 0.5" seeing. 

For the reasons mentioned in Sect. l3.2l we preferred to plot 
the modified signal to noise S/N* vs. My instead of the 
canonical S/N, as S/N* takes into account seeing effects 
and thus allows a comparison between measurements 
obtained with different seeing. A 2nd order polynomial 
fit to S/N* vs My is plotted over the data points. 
The corresponding fit for the canonical S/N is given as 
well. It is interesting to note that for 0.5" seeing the 
modified S/N* is consistently higher than the canonical 
S/N, while for 1.0" it is the opposite. This shows that 
for our simulations, 0.75" is about the limiting seeing 
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Fig. 2. Examples of simulated dEs in the / band. For the three upper examples, the image sequence is: Original image 

— model subtracted image — Power spectrum of model subtracted, square root divided and cleaned image. For the 
lowest example, the sequence is: Model subtracted image — Cleaned model subtracted and square root divided image 

— Power spectrum of the latter image. The different sequence for the lowest image was chosen to illustrate masking of 
contaminating point sources. For the power spectra, the dashed line is the best fit of P(k) — PSF(k) x Pq + P\, when 
rejecting the points marked with crosses. For technical reasons, the SBF amplitude at 0.5 wavenumber is equal to the 
amplitude at wavenumber 1.5. For the first example from the top, the dotted line represents a fit to all datapoints 
except the outermost and innermost one. The integration time was 1 hour for all examples. Parameters of the simulated 
dEs: 1st panel from top: My = —10.57 mag, (m — M) = 29.4 mag, /iy=25.5 mag arcsec^ 2 , (V — J)=0.88 mag, 
M\=-2A mag, dmi = mi, s imuiatcd — mi, measured = —0.22 mag. Seeing 0.5". 2nd panel from top: My — —11.13 mag, 
(m - M) = 31.4 mag, ^ v =25.1 mag arcsec~ 2 , (V - J)=0.90 mag, Mt=-2.86 mag, Srrii = 0.06 mag. Seeing^". 3rd 
panel from top: M v = -13.72 mag, (m — M) = 31.4 mag, ^=23.42 mag arcsec" 2 , (V - 7)=0.99 mag, Mt=-2.44 
mag, 8m\ = —0.02 mag. Seeing 1.0". 4th panel from top: My = —16.74 mag, (m — M) = 33.4 mag, fiy—20.9 mag 
arcsec~ 2 , (V - 7)=1.10 mag, Mt=-1.98 mag, Srrii = -0.03 mag. Seeing 0.5". 

below/above which S/N* becomes higher/lower than the canonical S/N. 
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Fig. 6. Results of the simulations for (m — M)=31.4 and seeing=0.5" at varying integration time / zero point 
ZP=27.0+2.5*log(t). From left to right: t= 900; 1800; 3600; 7200 seconds. 



One can see that for the highest S/N* data, the relative 
accuracy of the SBF method is of the order of 0.1-0.2 mag, 
i.e. 5-10% in distance. By glancing at the bottom panels 
of Figures |21 to El one immediately notices an average 
offset between simulated and measured SBF amplitude of 
the order of 0.15 mag in the sense of measuring too faint 
SBF. The mean offset, when disregarding the two most 
extreme offsets, is 0.14 mag. This offset is dealt with in 
more detail in Sect. IP 



4.1. Limiting absolute magnitudes 

For each of the 9 simulated sets of dEs a limiting absolute 
magnitude Mv.iim was determined, below which the deter- 
mination of cluster membership is not reliable anymore. 
These are the two conditions for reliable cluster member- 
ship determination we adopted: 

First, the difference between offset corrected measured 
and simulated mi must be smaller than 0.5 mag. The mean 
offset with regard to which the measurement difference is 
dchned is denoted as Sm\ and indicated in Figures |3| to (BJ 
The value of 0.5 mag was chosen as it is about equal to 
the maximum uncertainty in deriving M\ from (V—I), see 
Sect. |21 To require a higher measurement accuracy than 
the intrinsic uncertainty of the method would be unnec- 
essary. 

Second, the modified S/N* of the measurement must be 
higher than 6. This limit was adopted to avoid an SBF 
measurement of an object with S/N* < 4 mimicking a 
S/N* of the order of 5 or 6 because of a measured mi 0.5 
mag brighter than simulated. 

Mv,iim is then defined as the absolute magnitude at which 
50% of the measured galaxies fulfil the above criteria. It is 
indicated for each set in Figures|3to|ni As 5mi and Mv,ii ro 
depend on each other, they were determined iteratively. 





0.5" 


1.0" 


29.4 


-9.7 


-10.1 


31.4 


-12.8; -11.7; -10.8; -10.8 


-12.75 


33.4 


-14.9 


-16.75 



Table 1. Limiting absolute magnitude Mv,u m for the dif- 
ferent combinations of distance modulus (denoted in the 
left column) and seeing (denoted in the top line). The 
integration time t was 3600 seconds except for the four 
different values at 0.5" and 31.4 which correspond to four 
different integration times t=900, 1800, 3600 and 7200 sec- 
onds from left to right. The zero point including integra- 
tion time t is given by 27 + 2.5 * log(t). 



In Table ^ -^vjim is tabulated. One can see that for 
the two smaller distance moduli 29.4 and 31.4, the SBF- 
Mcthod can reach very faint magnitudes. For galaxies with 
My — —10 to —11 mag within a distance of about 20 
Mpc, reliable SBF measurements with accuracies better 
than 0.5 mag can be obtained at about 0.5" seeing, the 
given zeropoint of 27 mag and an integration time of lhr. 
The SBF-Method is therefore a very valuable tool for ex- 
tragalactic distance measurements, even for the faintest 
galaxies. 

What are the effects of varying seeing, integration time 
and distance modulus? 



4.1.1. Effects of varying seeing 

Table shows that for the distance moduli 31.4 and 
33.4, increasing the seeing by a factor of 2 brightens 
Mv,iim by about 2 mag. This corresponds to about 1.4 
mag in central surface brightness /Uo(V), or a factor of 
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Fig. 3. Results of the simulations for (m — M)=29.4 mag 
and zero point ZP = 27.0 + 2.5 * log{3600) mag. Left 
panel: Seeing=0.5". Right panel: Seeing=1.0". Top 
panel: Logarithm of S/N* (definition see text) plotted 
vs My. The dashed line is a 2nd order polynomial fit to 
the data points. The dotted line is the fit corresponding 
to the canonically defined S/N (see text). Bottom panel: 
8mi = TOLsimuiated - 7ni,mca S urcd plotted vs. My. For mag- 
nitudes fainter than the limiting magnitude My indicated 
by the dotted vertical line, more than 50% of the mea- 
sured galaxies have S/N* < 6 or a deviation of more than 
0.5 mag from the horizontal dashed line. The latter one 
denotes the mean Srfii for the measured galaxies brighter 
than My. 8mi and My were determined iteratively. 
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M v [mag] 
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Fig. 4. Results of the simulations for (m — M)=31.4 mag 
Symbols/lines as in Fig. EI 



Fig. 5. Results of the simulations for (m — M)=33.4 mag. 
Symbols/lines as in Fig. EI 



3.6 in central intensity. This increase by almost a factor 
of 4 in central intensity is plausible: the pixel-to-pixel 
fluctuations are smoothed by a factor of 2 at 1.0" seeing 
compared to 0.5" while the SBF amplitude is proportional 
to the square root of the intensity; therefore the intensity 
must be increased by a factor of 4 to compensate for the 
smoothing caused by the twice as large seeing. 
The "rule" extracted from that behaviour is: Increasing 
seeing by a factor x needs increase of intensity by a factor 
of x 2 to be compensated. 

For the 29.4 mag distance modulus, there is no big 
difference between 0.5" and 1.0" seeing. This is because 
at the faint magnitudes around —10 mag, the central 
surface brightness is about 25 mag arcsec -2 in /, which 
is only 2-3 er above the sky standard deviation, i.e. the 
mean surface brightness is close to the detection limit, 
and measuring its fluctuations is very difficult, even if 
the SBF amplitude is not much smaller than the mean 
surface brightness. That is why going from 1.0" to 0.5" 
seeing, no significant improvement of limiting absolute 
magnitude is reached for the 29.4 mag distance modulus. 



4.1.2. Effects of varying integration time / zero point 

What is the necessary scaling in integration time t 
to account for varying seeing? At 31.4 mag distance 
modulus, only <=900 seconds are needed for 0.5" seeing 
to reach the same Myjim as for 1.0" seeing and £=3600 
seconds. 

Thus, increasing seeing by a factor x needs increase of 
integration time by a factor of x 2 to be compensated. 
This is equivalent to keeping integration time fixed and 
increasing the zero point by 2.5*log(a; 2 ). 
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Tabic n shows that increasing the integration time t by 
a factor of 2 results in a 1 magnitude fainter Mv,Um) or 
about a factor of 2 in central intensity. As S /N and SBF 
amplitude are proportional to \fi and V7, respectively, 
this result should be expected. 

Thus, increasing integration time t by a factor x allows 
SBF measurement for objects with central intensity 
fainter by the same factor x. This is equivalent to keeping 
integration time fixed and increasing the zero point by 
2.5*log(a; 2 ). 

No notable change in limiting magnitude is seen when 
increasing t from 3600 to 7200 seconds. This might 
be partially due to statistical reasons, but the major 
reason is that the mean surface brightness is close to 
the detection limit (about 5 sigma above the sky noise) 
and the angular extent of the simulated dEs is only a 
few arcseconds. For galaxies with My ~ —11 mag, the 
additionally detected region when going from 3600 to 
7200 seconds carries no measurable SBF signal anymore. 



4.1.3. Effects of varying distance modulus 

The strength of the SBF relative to the underlying 
mean surface brightness decreases linearly with distance. 
As the SBF are proportional to the square root of the 
intensity, the intensity must increase by a factor of x 2 
when distance increases by a factor of x to compensate 
for that. Table shows that increasing distance modulus 
by 2 mag results in a 3-4 mag brighter limiting magnitude 
Mv.iim- This corresponds to about a factor of 10 in 
central intensity, which is slightly more than the expected 
value of 2.5 2 =6.25. The reason for this is that the angular 
area over which the SBF signal is sampled is smaller at 
2.5 higher distance for the same object. 



4.1.4. A simple rule 

Summarizing the scaling relations found in the last three 
subsections, we give the following rule to calculate the 
limiting magnitude M Vnew at new observing conditions 
different to the reference ones adopted in our simulations: 

M£ new = M£ rof -6.64*Zo5(— )-1.33*(ZP rcf -ZP now )(7) 

-^ref 

My Io{ is the limiting magnitude calculated from our 
simulations at the given distance modulus. s re f is the 
seeing diameter in our simulations, s new the new seeing 
diameter. ZP Te { is the total zero point in the simulations, 
i.e. ZP rci = 27.0 + 2.5 * log(t) with t being the total 
integration time. ZP ncw is then the new total zero 
point. This all refers to a gain of 1, i.e. the zero point is 
expressend in terms of electrons and not ADU. 
Note, however, that equation[7|is restricted to cases where 
the mean surface brightness of the galaxy is significantly 
higher than the sky noise. As when the surface brightness 



gets too close to the sky noise (less than about 5 sigma, see 
the former three subsections), changing integration time 
or seeing does not have strong effects on SBF detectability. 



4.2. A bias in the simulations 

To find out the reasons for the 0.14 mag average offset 
between simulated and measured mi, a number of tests 
were performed: 

First, an area of constant surface brightness and SBF 
amplitude was simulated, but not convolved with the 
seeing. After subtracting the mean brightness, dividing 
by its square root and calculating the power spectrum, 
the resulting image should have a mean equal to the SBF 
amplitude. This was the case to within 1%, independent 
of surface brightness and SBF amplitude. 
Second, an area of intensity zero was simulated with only 
one pixel given an intensity Iq different from zero. This 
area was then convolved with the seeing profile. The mean 
of the resulting image was 0.99 x Iq. he. by restricting 
the seeing psf simulation to 7 times the FWHM, about 
1% flux is lost. This is because we chose a realistic Moffat 
profile, which has stronger wings than the commonly 
used Gaussian. 

The two above tests show that the algorithm used 
to implement and measure pixel-to-pixel SBF works 
correctly, and that convolution with the seeing profile to 
7 times the FWHM imposes a negligibly small flux loss. 
The final check performed was to simulate areas of several 
hundred pixel side length with constant surface brightness 
and SBF amplitude, which were then convolved with the 
seeing. After subtraction of the mean brightness, division 
by its square root, calculation of the power spectrum and 
azimuthally averaging, the SBF amplitude according to 
equation (@J - i.e. in the limit of k=0 - was too faint 
by 0.10 to 0.15 mag, independent of the strength of the 
simulated SBF, and for both seeing values 0.5 and 1.0". 
This offset corresponds to the one found in the main 
simulations. It shows that recovering the underlying 
pixel-to-pixel fluctuations from seeing convolved images 
is subject to small, but non-negligible loss of fluctuation 
signal, at least in our simulations. 

Already Tonry & Schneider ( 1988J noted a bias towards 
measuring too faint SBF of the order of 10% in distance 
or 0.2 mag in distance modulus, and attributed this to 
truncation of the seeing profile in their simulations. As 
mentioned above, in our case seeing truncation imposes a 
negligible flux loss. 

Zero point problems of the order of 0.15 mag are certainly 
a matter of concern when one aims at determining 
absolute values like Hq. However, our aim was to deter- 
mine limiting magnitudes and surface brightnesses for 
cluster membership determination of dEs when using the 
SBF-Method. These measurements all rely on relative 
distances and are therefore independent of any zero 
point offset. The offset found in our simulations does not 
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Fig. 7. Apparent SBF magnitude rn\ plotted vs. (V — I)q 
for 6 Centaurus cluster dEs (from Mieske et al. |in prep.| ) . 
The dashed line is a fit to the data when adopting a slope 
of 4.5, as in Tonry's equation The corresponding dis- 
tance modulus is 33.15 ± 0.04 mag. When fitting slope and 
zero point at the same time, the resulting slope differs by 
less than one a from 4.5. 

significantly change the statements made in the former 
sections. 

We note that unresolved background galaxies generally 
increase the measured SBF signal. However, for our simu- 
lations this contribution is negligible. Using formula (13) 
of Jensen et al. (1998) for the /-band and inserting the 
values used for the magnitude distribution of background 
objects, we get for the relative contribution of background 
galaxies to the SBF signal at 33.4 distance modulus a 
value of the order of 2-4%, depending on seeing and the 
galaxy's magnitude. For the distance moduli 29.4 and 
31.4, the contribution is below 1%. 




15 15.5 16 16.5 17 17.5 
V o [ ma g] 

Fig. 8. S/N of the SBF measurement vs. apparent Vq. 
Large circles are the real Centaurus data (from Mieske 
et al. |in prep] ). Dots represent galaxies simulated accord- 
ing to the Centaurus data's colour-SBF relation, colour- 
magnitude relation and surface brightness-magnitude re- 
lation. The dashed line represents a fit to the real data, 
the solid line is a fit to the simulated data. Both fits are 
consistent with each other within the error ranges. 

real measurements are plotted over the results for the sim- 
ulations. A line is fit to both real and simulated data. One 
can see that the simulations do not overestimate the S /N 
of the real data. Both fits are consistent with each other. 
This consistency between real and simulated data confirms 
that the simulations presented in the previous sections are 
a good approximation of reality. The applied generaliza- 
tions like zero ellipticity and purely exponential profile 
apparently do not introduce a notable bias towards too 
high or too low S/N. 



4.3. Comparing real SBF data with simulations 

The S/N achieved in SBF measurements for 6 bright 
Centaurus cluster dEs from VLT FORS1 images is com- 
pared with the S/N obtained from simulations tuned to 
reproduce the measured values. 

Fig.[7|shows the colour-SBF relation found for the 6 galax- 
ies (Mieske et al., |in prep.| ). They cover a magnitude range 
of 15.4 < Vo < 17 A mag. The colour-SBF relation from 
Fig. 13 a colour-magnitude relation and surface brightness- 
magnitude relation were fit to the measured values of the 
Centaurus dEs, and 64 galaxies were simulated according 
to these relations. Their SBF amplitude was measured as 
described in Sect. El In Fig. OH the log(S/N) values of the 



5. Summary and conclusions 

Extensive simulations of SBF measurements on dEs for 
three different distance moduli 29.4, 31.4 and 33.4 mag, 
two different seeings 0.5" and 1.0" and four different 
observing times 900, 1800, 3600 and 7200 seconds have 
been presented. For each of the simulated sets of dEs, the 
limiting magnitude My.ii m below which a distance mea- 
surement is not reliable anymore has been determined. It 
was shown that for distances < 20 Mpc, the SBF method 
can yield reliable cluster membership of dEs down to 
very faint limiting magnitudes, e.g. Myjim — —10 mag 
for a distance of 7.5 Mpc, and Afv,iim — —11 mag for 
19 Mpc distance, at lhr integration time, 0.5" seeing 
and a zero point of 27 mag in the /-band. For the SBF 



S. Mieske et al.: SBF distances to dEs in nearby clusters 



11 



measurements, a modified signal to noise S/N* has been 
defined that incorporates the seeing dependence of SBF 
detectability. 

The effects of varying seeing, integration time and 
distance modulus on the limiting absolute magnitude 
are investigated. A number of simple rules, confirmed 
by theoretical considerations, are derived in order to 
calculate limiting magnitudes, needed integration times 
or seeing for observing conditions different to the ones 
adopted for our simulations. 

It is pointed out that the total uncertainty in obtaining 
a distance modulus for a dE with SBF measurements is 
the quadratic sum of the measurement uncertainties for 
mi and the uncertainties in deriving Mi from (V — I). 
As both uncertainties are 0.5 mag or smaller, the worst 
possible measurement accuracy is of the order of 0.65 
mag or 35% in distance. This would apply to very 
faint and blue dEs close to the limiting magnitude. For 
brighter and redder dEs, the total uncertainty decreases 
significantly, as deriving Mi from (V — I) is less uncertain 
for red objects and the measurement accuracy improves 
to about 0.1 mag for the brightest simulated dEs. While 
the uncertainty in deriving Mi from (V — I) imposes a 
lower limit on the distance accuracy when observing field 
dEs, it allows rough age-metallicity estimates for blue 
cluster dEs, as most of the times the cluster is separated 
from background/foreground galaxies by significantly 
more than 1 mag in distance modulus. 
We find that on average the measured SBF magnitude 
is 0.15 mag fainter than the simulated one. A number of 
tests show that this is due to loss of fluctuation signal 
when recovering pixel-to-pixel fluctuations from a seeing 
convolved image. This average offset shows that much 
care must be taken when deriving absolute values like Hq 
with the SBF-Mcthod. When aiming at relative distances 
like for cluster membership determination, a bias is of no 
concern. 

By comparing real SBF data of Centaurus Cluster dEs 
with simulations tuned to reproduce the real data, we find 
that our simulations do not overestimate the achievable 
S/N of the SBF method, but are consistent with real 
measurements. Therefore the statements about limiting 
magnitudes for the technique made in this paper are 
reasonable and we would not expect a very different 
behaviour in real observations. 

An ideal application of the SBF technique would be a 
deep and wide field survey of several nearby clusters such 
as Fornax, Virgo or Doradus. With the arrival of wide 
field cameras on large telescopes (Suprime cam on the 
Subaru telescope or IMACS on Magellan), this is a very 
promising possibility to determine well the very faint end 
of the galaxy luminosity function in nearby clusters. 
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